Blazar halos as probe for extragalactic magnetic fields and maximal acceleration 
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High energy photons from blazars interact within tens of kpc with the extragalactic photon 
background, initiating electromagnetic pair cascades. The charged component of such cascades is 
deflected by extragalactic magnetic fields (EGMF), leading to halos even around initially point- 
like sources. We calculate the intensity profile of the resulting secondary high-energy photons for 
different assumptions on the initial source spectrum and the strength of the EGMF, employing also 
fields found earlier in a constrained simulation of structure formation including MHD processes. We 
find that the observation of halos around blazars like Mrk 180 probes an interesting range of EGMF 
strengths and acceleration models: In particular, blazar halos test if the photon energy spectrum at 
the source extends beyond ~ 100 TeV and how anisotropic this high energy component is emitted. 

PACS numbers: 98.70.Sa, 95.85.Pw, 98.54.Cm, 98.62. En 
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I. INTRODUCTION 

Observations of photons with energies up to tens 
of TeV, mainly by imaging air Cherenkov telescopes 
(I ACT), have opened in the last decade a new window 
on the non-thermal side of the Universe [l[ . One of the 
main aims of these observations is the identification of 
the sources of cosmic rays and their acceleration mecha- 
nism. Observations of very high energy photons can how- 
ever provide information on many more areas of astro- 
physics, as e.g. the amount of extragalactic background 
light (EBL) and the strength of extragalactic magnetic 
fields (EGMF). 

The tight relationship between the EBL, the strength 
of EGMFs and the horizon of high energy photons is 
well-known: Photons with energy above the pair creation 
threshold, E th ~ 30 TeV for E ~ 10~ 2 cV as characteris- 
tic energy of the background photons, can interact with 
them, initiating electromagnetic pair cascades 0- The 
charged component of such cascades is in turn deflected 
by magnetic fields, first by fields in the cluster surround- 
ing the source and then by EGMFs Q. If the EGMF 
is not extremely week, B < 10~ 16 G, secondary photons 
are deflected outside a point-like source and thus the ob- 
served point-like photon flux 1(E) consists only of the 
surviving primary photons, 1(E) = exp[— t 7J (E)]Iq, with 
r 77 as the depth for pair production. If the EGMF ex- 
ceeds 1CU 10 G in a large volume fraction, the halo be- 
comes too large and its intensity therefore too small for 
a detection. In the intermediate range, (10~ 16 -10~ 10 ) G, 
the deflections and thus the halo extension may be suf- 
ficiently small that an observation by current or next 
generation IACTs may be feasible Q ■ 

We focus in this work on the effect of EGMFs on elec- 
tromagnetic cascades, providing detailed predictions for 
the expected intensity I(E, d) of the halo produced by 
cascades as function of the observed energy E and the 



angular distance i) to the source. Since the EGMF is 
highly structured and the Milky Way is contained in the 
supergalactic plane, the magnetic field along the line-of- 
sight to different blazars, and thus also their halos, will 
vary significantly, even if they would be at the same dis- 
tance and would have the same source energy spectrum. 
It is therefore crucial to use an EGMF model that de- 
scribes at least qualitative correctly the location of voids 
and filaments. Since deflections close to the observer are 
most important, it is sufficient to account for the EGMF 
structure of the local universe. For our predictions of the 
expected halo intensity I(E, fi) around selected blazars 
we use therefore EGMF models calculated in Ref. [f| that 
are based on a constrained simulation of the large-scale 
structure within w 115Mpc. 

Our results confirm the estimate of Ref. [1] that blazar 
halos are detectable for specific acceleration models, if 
EGMFs are as weak as found in the simulations of 
Ref. [f|. Thus blazar halos may provide valuable infor- 
mation about the strength of EGMFs in a range where 
other tools like Faraday rotation measurements can give 
only upper limits. Moreover, our detailed predictions 
for the expected intensity I(E, t?) of the halo produced 
by the electromagnetic cascade allow one to disentangle 
the contributions from secondaries generated close to the 
source, as suggested in Ref. and during propagation 
in the EGMF. The intensity of the halo is also sensitive 
to the total luminosity of the source above ~ 100 TeV 
and to the opening angle ?9 of the cone in which most of 
this high energy radiation is emitted. 

We find that the measured energy spectrum of most 
blazars detected with IACTs may be described by a 
broken power-law that is steeper in the measured en- 
ergy range (lOOGeV-few TeV) than usually assumed, 
if one adds a softer high energy component. This ad- 
ditional high-energy component is absent both in Syn- 
chrotron Sclf-Compton SSC and Inverse Compton mod- 
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els for blazars [|| , but may be caused by hadrons @, § • 
In most hadronic models, e.g. acceleration close to the 
core in electromagnetic fields or in hot spots, the bulk of 
radiation, although emitted anisotropically, extends over 
a cone with opening angles aj et > 10°. For such condi- 
tions, we found that the halo is not strongly suppressed. 
Thus the observation of blazar halos would be an indi- 
cation for an extension of the source photon spectrum 
beyond 100 TeV that is not strongly beamed and thus 
also evidence for the acceleration of hadrons in blazars. 

This article is structured as follows: We start in Sec.HTl 
with a brief description of our Monte Carlo (MC) scheme. 
Then we study in Sec. IIIII a toy model that allows us to 
exhibit clearly the dependence of the halo properties on 
external parameters as the strength of the EGMF, the 
source photon energy spectrum and the amount of beam- 
ing. In Sec. lIVl we use magnetic fields from Ref. [f| along 
the line-of-sight towards four selected blazars and cal- 
culate their halos, choosing the source energy spectrum 
such that the observed photon spectra at low energies are 
reproduced. Finally, we summarize in Sec. |Vj 

II. MONTE CARLO SCHEME 

We choose to describe the propagation of the electro- 
magnetic (e/m) cascade in the extragalactic space using a 
MC method. Starting from an initial photon emitted by 
the source, we trace the development of the cascade tak- 
ing into account e + e~ pair production by photons and in- 
verse Compton scattering (ICS) of electrons due to their 
interactions with photons from the radio, infra-red, and 
the microwave background as well as synchrotron energy 
losses of electrons. For the spectral densities of the back- 
grounds we use the "Best-Fit06" model from Ref. @ for 
the EBL and Ref. for the radio background. 

Since we are interested in gamma-rays arriving to the 
observer with small angular deflections, we can use an 
essentially one-dimensional picture for the propagation 
of the e/m cascade. A typical propagation pattern is 
sketched in Fig. [TJ The initial photon is emitted by the 
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FIG. 1: Schematic view of the electromagnetic cascade in 
extragalactic space. 

source 5" with the angle a with respect to the direction 
to the observer O. Due to interactions of both the initial 
photon and of secondary leptons, an e/m cascade devel- 
ops over the distance x, until at point P the final photon 



is produced. Note that in the ultrarclativistic limit rel- 
evant for us, the daughter particles move in the forward 
direction. The final photon arrives to the observer under 
the angle $ with respect to the direction to the source, 
while its wave-vector forms the angle (3 with respect to 
the initial photon direction. The three angles i9, a, and 
(3 are related by 

a = p-fi, (1) 

X 

sin$ = — sin/? . (2) 

For small $ we are interested in, we have 

sin/3. (3) 

We calculate the average scattering angle f3 within the 
small angle approximation as 



with (3f as angular deflection of the electron i during the 
cascade development. The deflection angle in a regular 
magnetic field is 




for a particle with momentum p± perpendicular to B 
propagating the distance d. 

As it can be read off from Eq. ([3|) , the effect of deflec- 
tions (3 at large distances from the observer is suppressed 
by the ratio x/L compared to those close by. Hence an 
important contribution to the halo at small $ are final 
photons produced at x <C L, and passing the remained 
distance L — x without interactions. A competing al- 
though less important contribution comes from photons 
with small scattering angle (3. In particular, the latter is 
somewhat enhanced in the case of a structured EGMF: 
Filaments of the EGMF may be crossed by intermediate 
photons in the e/m cascade, while intermediate electrons 
propagate in the voids. As a result, the strength and the 
structure of the magnetic field close to the observer is 
more important than the "average" field. 



III. PARAMETER DEPENDENCE OF THE 
HALO IN A TOY MODEL 

The intensity and the extension of the halo around 
a blazar caused by the deflection of e + e~ pairs in the 
electromagnetic cascade depends mainly on the strength 
and the structure of the EGMF, the spectral shape of the 
source spectrum and the total luminosity of the source 
above ~ 30 TeV. In order to disentangle better the influ- 
ence of these different factors, we employ in this section 
a toy model: We use a uniform extragalactic magnetic 
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FIG. 2: The differential energy spectrum I(E, ■&) for three dif- 
ferent angular bins and a uniform magnetic field with strength 
B = 10" 12 G (black), 10" 11 G (blue) and B = 10" 10 G (red). 



field and a simple power-law as the energy spectrum of 
photons leaving a blazar. 

If not otherwise specified, we use as reference param- 
eters a uniform EGMF with field strength B = 10" 12 G, 
a power-law for the source energy spectrum of photons, 
dN/dE oc E'" 1 , with 7 = 2, £ max = 10 17 cV and place 
our reference source at the distance of 135 Mpc, i.e. rep- 
resentative for Mrk 501. 

a. Dependence on the EGMF strength We first dis- 
cuss the change of the intensity profile I(E, fl) = 
d 2 N(E,ti)/drt varying the strength of the EGMF. Fig- 
ure [2] shows the differential energy spectrum I(E, •&) as 
function of the observed energy E for three i? bins: The 
solid lines include both the point-like and the halo com- 
ponent within d < i?o = 0.06° (from now on called 'point- 
like'), while the dotted lines denote the mean intensity 
between do < "& < 2$o ('halo 1') and the dashed lines de- 
note the intensity at $ = 1.2° ('halo 2'). Our choice for 
the angular binning is related to the angular accuracy 
of present IACT (typically, few arcmin) and to experi- 
mental methods of studying blazar halos: The intensity 
of 7-radiation from the angular vicinity of the source is 
compared to the background intensity in the field of view 
of the apparatus (typically within 1 degree). Thus, the 
feasibility of halo measurements can be characterized by 
the relative difference in the integral intensity between 
'halo 1' and 'halo 2'. On the other hand, the overall sen- 
sitivity of an instrument to the halo component can be 
determined from the intensity ratio between the 'point- 



like' and 'halo 1' components. As the strength of the 
EGMF increases from B = 10" 12 G (left) to B = 10" 10 G 
(right), the /(> Eq,{>) distributions shown in Fig. ^flat- 
ten, since the deflection angles increase. Moreover, the 
/(> Eo,-d) distributions for different Eo become more 
similar to each other. As a result, the ratio between 
the halo intensity at small ■& and the point-like intensity 
decreases too. As one can see from the figures, even for 
relatively strong magnetic fields with B ~ 10" 10 G a halo 
observation is possible, in principle, if one restricts one- 
self with gamma energies in excess of 10 TeV. However, 
imposing such a cutoff would result in a drastic reduction 
of the observed flux, the measurement thus being hardly 
possible with present IACT sensitivity. Thus, the obser- 
vation of blazar halos is more likely, if the line-of-sight 
to the blazar is contained in a void, where the fields of 
strength B = 10" 12 G may be typical. 

b. Dependence on the spectral slope A comparison 
of the intensity profiles for two different values of the 
exponent, 7 = 2.0 and 2.4, in the power-law for the 
source energy spectrum of photons, dN/dE ex _E" 7 , is 
shown in Fig. 3] The halo component depends more 
weakly on the spectral slope 7 than the point-like com- 
ponent, but its relative normalization decreases sharply 
for steeper source energy spectra. As a result, the ratio 
of halo and point flux decreases strongly towards lower 
energies. Nevertheless, the total number of halo events 
increases lowering the threshold energy Eq. The opti- 
mal energy to look for the halo component depends also 
on the background of cosmic ray events and the back- 
ground of point-like events with misreconstructed arrival 
direction. Therefore detailed studies taking into account 
the specific detector properties are required for such an 
optimization. 

c. Dependence on the distance We compare the dif- 
ferential energy spectrum I(E, $) of the point-like to the 
'halo 1' and 'halo 2' components for two different dis- 
tances in Fig. If one factors out the trivial distance 
dependence, i.e. normalizes the flux at say 1 TcV to a 
constant value, then the halo components of sources at 
135 and 300 Mpc show only minor differences: This be- 
havior is explained by the fact that deflections at the 
distance I to the observer are weighted with the factor 
~ (L — l)/L. Hence, deflections close to the observer 
contribute most to the size of the halo, while interactions 
at large distances lead mainly to a reduction of the ob- 
served flux. Moreover, the strongest deflection in case of 
a uniform EGMF considered here comes from the propa- 
gation of the last electron in the e/m cascade, which has 
the smallest energy, hence, is most strongly deflected by 
the field. 

d. Dependence on high-energy luminosity We use 
again a power-law dN/dE oc E~ 2 , varying now the max- 
imal energy £ max from £? max = 10 13 , 3 x 10 13 , 10 14 to 
10 15 eV. In Fig.[6l the influence of the luminosity injected 
above the pair-production threshold on the halo intensity 
is clearly visible. While close to the threshold an increase 
of -E max leads to a strong enhancement of the halo com- 
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FIG. 3: The integral intensity J(> Eo, ■&) as function of the angular distance $ to the source for a uniform magnetic field with 
strength B = 10" 12 G (left) and B = 10" 10 G (right). 




FIG. 4: Comparison of ^-dependence of integral intensity /(> Eo,$) (left panel) and of differential energy spectra I(E, , &) 
(right panel) for different source spectra: power-laws with 7 = 2.0 and 7 = 2.4 (lower curves in the right panel; solid lines for 
point-like, dotted for 'halo 1', and dashed for 'halo 2'). A uniform magnetic field with strength B = 10 -11 G and the distance 
134 Mpc is assumed. 



ponent, a further increase of E max beyond 10 15 eV has 
much smaller influence of the halo intensity for the cho- 
sen 1/E 2 spectrum. 

e. Dependence on the jet opening angle ev max Fi- 
nally we discuss the influence on the halo intensity of an 
anisotropic emission of the high energy photons. Qual- 



itatively, one expects from Fig. 1, that restricting a to 
smaller and smaller values for a given 1} requires that the 
last interaction at P happens closer and closer to the ob- 
server. Thus the intensity of the halo will be reduced on 
all angular scales, but we do not expect an angular cutoff 
as effect of an anisotropic emission by the source. 
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FIG. 5: Comparison of differential energy spectra for three 
angular bins for 135 and 300 Mpc source distance. 



In order to obtain a quantitative understanding for 
the importance of beamed emission, we have recalculated 
our standard case (isotropic gamma-ray emission in the 
source) restricting possible emission angles a with respect 
to the line-of-sight to the blazar to Oj C t = 60°, 30°, 15° 
and 5°. The resulting integral intensity > 300GeV) 
as function of $ is shown in the left panel of Fig. while 
we show the differential energy spectrum I{E, •&) of the 
point-like to the 'halo 1' and 'halo 2' components in the 
right panel. The integral intensity > 300GcV) shows 
clearly an overall suppression as function of i? instead of a 
cutoff. The extent of this suppression becomes smaller at 
higher energies, but clearly the detection of blazar halos 
would disfavour the emission of the high energy compo- 
nent in a very narrow cone, aj e t ^ 10°. 

In summary, we have found that the shape of the halo 
depends sensitively on all investigated parameters, but 
generally I(E, , d) decreases faster than the l/i? behavior 
suggested in Ref. Q ■ The obtained enhancement of the 
small d region, compared to Ref. [H , is due to the precise 
MC description of the cascade process, being related to 
the contributions from the tails of the distribution of the 
electron's interaction length and of the secondary photon 
energies in ICS: Shorter distances travelled by electrons 
and higher energy of the parent electron producing the 
final photon result in smaller than average values of 
Similarly, decreasing the opening angle of the high energy 
emission of photons in the source does not produce an 
angular cutoff in the observed photon spectrum, contrary 
to what has been suggested by the authors of Ref. [J]. 
The latter was an artifact of considering only the first 
two interaction steps in the electromagnetic cascade in 



FIG. 6: Comparison of the differential energy spectra for the 
point-like (three solid lines at the top) and "halo 1" (dashed 
lines) components, with three different energy cutoffs of the 
source spectrum: top down -E max = 10 15 ,10 14 ,3 x fO 13 and 
10 13 eV. 



Ref. [4(. Most importantly, we found that the intensity 
of the halo can exceed 1% of the point-like contribution 
at energies in excess of few TcV for field strengths B = 
l(r 12 G. 



IV. HALO PREDICTIONS FOR SELECTED 
BLAZARS 

We now turn from the toy model studied in the last 
section to a realistic study of the blazar halos for the five 
nearest blazars. The table |T] summarizes the distance, 
the observed (7) and intrinsic (7^) spectral index assum- 
ing a power-law, the energy of the observed cut-off or 
the energy of the last measured data point as well as 
the reference used. As far as possible we used the pho- 
ton spectra measured during quiet periods. While in the 
case of Mrk 421 the maximal energy E max = 1.4 TeV 
corresponds to an observed cutoff in the measured en- 
ergy spectrum, the numerical value of E max corresponds 
otherwise to the last data point given in the cited refer- 
ences. From our discussion of the toy model, it is clear 
that the low luminosity of Mrk 421 above the pair cre- 
ation threshold makes it an unfavorable candidate for the 
search of blazar halo. Therefore we will restrict our anal- 
ysis to the remaining four blazars, Mrk 180, Mrk 501, 
1ES 1959+650 and 1ES 2344+514. 

We used the EGMF along the line-of-sight towards the 
four blazars that results from a constrained realization 
of the local universe (see Ref. Q and references therein). 
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FIG. 7: Comparison of the ^-dependence of the integral intensity 7(> Eo, $) (left panel) and of the differential energy spectra 
I(E,$) (right panel) for four different jet opening angles, ajet = 60°, 30°, 15° and 5° (top down). 



In short, the initial density fluctuations were constructed 
from the IRAS 1.2-,Iy galaxy survey by smoothing the ob- 
served galaxy density field on a scale of 7 Mpc, evolving 
it linearly back in time, and then using it as a Gaussian 
constraint (see [l9j ) for an otherwise random realization 
of the ACDM cosmology. The volume constrained by the 
IRAS observations covers a sphere of radius w 115 Mpc 
centered on the Milky Way and therefore encloses the 
prominent local structure including a fair representation 
of the large scale density field. Many of the most promi- 
nent clusters observed locally can therefore be identified 
directly with halos in the simulation, and their positions 
and masses agree well with their simulated counterparts 
(see also [2(|). Although the smoothing of the observed 
galaxy density field when producing the initial condi- 
tions invokes spacial uncertainties within the simulated 
structures up to the used smoothing scale of 7 Mpc, the 
statistics of crossing filaments and voids within the local 
universe should nevertheless be quite realistic. This is es- 
pecially the case for the objects considered in this study 
as we can make use of the full path length through the 
simulated volume as these objects are lying outside the 
simulated region. 

For this work we used the magnetic field configuration 
obtained from two realizations (MHDy and MHDz) start- 
ing from different initial seed fields (for more details see 
Ref. H) having lower (MHDz) and higher (MHDy) ini- 
tial values for the magnetic seed field. Both simulations 
lead to magnetic fields within the galaxy clusters which 
are statistically in agreement with the still rare Faraday 
rotation measures and arc roughly bracketing the allowed 
range of magnetic fields in galaxy clusters. 



A large fraction of the line-of-sight towards the blazers 
of interest are within the volume covered by this sim- 
ulation, however it can cover only up to distances of 
ps 115 Mpc and we are forced to extend the magnetic field 
profile beyond that point. Since it is most important to 
reproduce the local field structure, we expect that the de- 
tails of the chosen extension have only a minor influence 
as long the extended field has the correct statistical prop- 
erties. In practise, we mirrored the field at the point of 
maximal distance along the line-of-sight. A comparison 
of the EGMF component perpendicular to the line-of- 
sight towards Mrk 180 obtained in the models is shown 
in Fig. [TO] 

We will first present results for the model MHDz that 
generally predicts stronger localisation of magnetic fields 
in cluster and filaments and thus smaller deflections than 
model MHDy. Later, we compare the differences between 
our results for the two models. 

We show our results 1 in Fig. [5] for Mrk 180 (top) and 
Mrk 501 (bottom) and in Fig. [U for 1ES 1959+650 (top) 
and 1ES 2344+514 (bottom). The differential energy 
spectrum of each source is calculated for two different 
injection spectra: In the first case, we assume that the 
observed spectrum extends with the same slope 71 up 
to E max = 10 17 cV. In the second case, we add a new, 
harder component with 72 < 71 that may be generated 
in particular by hadrons. Thus, we use as injection spec- 



1 Numerical tables of these results can be obtained from MK. 
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FIG. 8: Top: Mrk 180 (left panel: ^-dependence of integral intensity for the "broken" power-law initial spectrum with 71 = 3.5, 
7 = 1.8 and Eb = 2 TeV; right panel: differential energy spectra for three angular bins for the power-law source spectrum with 
7 = 2.8 and for the "broken" power- law case), bottom: Mrk 501 (left panel: ^-dependence of integral intensity for the "broken" 
power-law initial spectrum with 71 = 2.4, 7 = 2.0 and Eb = 5 TeV; right panel: differential energy spectra for three angular 
bins for the power-law source spectrum with 7 = 2.4 and for the "broken" power-law case); for all cases i?max = 10 17 eV. 



trum a broken power-law, adjusting the parameters like 
the break energy Eb such that we reproduce the observed 
spectral shape while optimizing the halo flux. 

A characteristic feature of the second case is a small 
bump between the break energy and the exponential cut- 
off at E t h introduced by the EBL interactions. Clearly, 
blazars which show a cutoff at energies below -Eebl can 
have only a minor additional hard component compared 
to blazars which show either no cutoff or even a turn up. 



As one can see from the Figures, for the "broken" power- 
law case we get very pronounced halos in the TeV range 
for Mrk 180 and 1ES 2344+514, with the halo intensity 
being comparable to the one of the point-like component. 

In Fig.QdJ we compare the energy spectrum of Mrk 180 
(left) and 1ES 2344+514 (right) obtained using the 
EGMF models MHDy and MHDz together with the "bro- 
ken" power-law initial spectra. For both sources, we ob- 
tain a reduction of both point-like and halo intensities for 
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FIG. 9: top: 1ES1959+650, (left panel: ^-dependence of integral intensity for the "broken" power-law initial spectrum with 
71 = 2.6, 7 = 1.8 and Eb — 3TeV; right panel: differential energy spectra for three angular bins for the power-law source 
spectrum with 7 = 2.58 and for the "broken" power-law case), bottom: 1ES2344+514 (left panel: ^-dependence of integral 
intensity for the "broken" power-law initial spectrum with 71 = 2.7, 7 = 1.7 and Eb = 5TeV; right panel: differential energy 
spectra for three angular bins for the power-law source spectrum with 7 = 2.7 and for the "broken" power-law case). 



the MHDy realization of EGMF which is characterized 
by stronger and less localized magnetic fields. As one can 
see from Fig. [TU1 the strength of the perpendicular com- 
ponent of EGMF in the MHDy realization is almost an 
order of magnitude higher than in the MHDz case. This 
results in a reduction of the halo component for both 
sources considered, especially at lower energies. Never- 
theless, for both realizations of the magnetic field struc- 
ture, Mrk 180 and 1ES 2344+514 remain very promising 



cases for the experimental detection of the halo in the 
TeV energy range. 



Finally we note that the HEGRA experiment derived 
an upper limit of 1% on the halo flux from Mrk 521 at 
0.5° during a burst period [l8[. 
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FIG. 10: The magnetic field B± perpendicular to the line- 
of-sight towards Mrk 180. 
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TABLE I: List of considered blazars. 



V. DISCUSSION 

We have studied blazar halos generated by deflections 
of the charged component of electromagnetic pair cas- 
cades. Our analysis has accounted both for the kinemat- 
ics of the pair-production and Compton processes with- 
out approximations as well as for the local structure of 
the EGMF by using EGMF models found earlier in a 
constrained simulation of structure formation including 
MHD processes. The intensity of these blazar halos de- 
pends mainly on the strength and the structure of the 
EGMF: Halos with detectable intensity require very weak 
EGMFs, comparable to those found in the MHDz simula- 
tions of Ref. Q . The direction of the EGMF vector close 
to the observer constitutes an uncertainty that cannot be 
resolved by numerical simulations. If the linc-of-sight of 
a specific source and the EGMF vector are aligned close 
to the observer, the effect of EGMFs will be reduced. 

Another necessary condition for the observability of 
halos is the extension of the source energy spectrum 
above ~ 10 14 eV. Such a high energy extension of the 
source spectrum may be generated by high energy pro- 
tons interacting with background photons in the source 
or emitting synchrotron or curvature radiation. Note 
however that the acceleration of hadrons in the source 



does not necessarily lead to a pronounced ultrahigh en- 
ergy tail of the accompanying gamma-ray spectrum. For 
instance, in the case of hadrons accelerated close to the 
core of the AGN, high energy gamma-rays start pair cas- 
cades on the background of ultraviolet photons. As a re- 
sult, the spectrum of gamma-rays leaving the source is es- 
sentially cut-off at 10 13 eV, for an example see e.g. Ref. [![. 
Thus, the maximal energies assumed in the present work 
for the photons leaving the source rather correspond to 
acceleration in larger volumes, filled with less dense pho- 
ton backgrounds and weaker magnetic field than close to 
the AGN core. Typical examples for such a scenario are 
the acceleration of protons in the outer or radio jets of 
AGNs. 

Throughout the paper, we implicitly assumed that the 
high energy part of the source spectrum has a hadronic 
origin and, consequently, is not strongly beamed. Since 
we found that the halo intensity is strongly suppressed 
if this high energy component is emitted too in a very 
narrow cone, aj e t < 10°, an observation of blazar ha- 
los would indicate a different emission mechanism for 
the source energy spectrum above ~ 10 eV, pointing 
to a hadronic origin. The obtained integral intensity 
I("d, > 300GcV) in the case of a beamed UHE emission 
shows clearly an overall suppression as function of the 
observation angle, without angular cutoff. The extent of 
this suppression becomes smaller at higher energies, but 
clearly the detection of blazar halos would disfavour the 
beamed character of high energy photon emission. 

The observed spectra of TeV blazar show strong time- 
variability. In contrast, any intrinsic time- variability of 
the halo component is washed out by time-delays in- 
duced by the curved trajectories of the charged secon- 
daries in the EGMF. Since the halo component has an 
hadronic origin in our scenario, we do not expect even an 
intrinsic time-variability to start with. Consequently, we 
used spectral data measured during 'quiet' periods for the 
sources considered. Similarly, experiments should com- 
pare the relative normalization of the halo to the 'point- 
like' component during quiet periods. 

In the most optimistic cases, the halo component found 
reaches 10% of the point-like component and may be al- 
ready detected or excluded by current IACTs. An obser- 
vation of blazar halos would provide valuable information 
on the EGMF along the line-of-sight to the source and 
on the minimal luminosity emitted above the pair cre- 
ation threshold. Moreover, a halo detection would point 
to an unbeamed emission of the high-energy part of the 
source photon spectrum, providing evidence for a differ- 
ent acceleration and emission mechanism than the stan- 
dard leptonic one. 

Finally, we comment on the generation of halos not by 
the electromagnetic cascade in the extragalactic space 
but close to the source. A calculation of the halo inten- 
sity from this alternative mechanism requires a detailed 
model for the photon backgrounds and the magnetic field 
in the source. Moreover, the model for acceleration, in 
particular the acceleration site, has to be specified. While 
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FIG. 11: Comparison of differential energy spectra for three angular bins for two realizations of EGMF (MHDy [red] and 
MHDz [blue]) for Mrk 180 (left panel) and for 1ES 2344+514 (right panel). 



it is thus clear that an additional contribution to the halo 
from processes close to the source exists, detailed predic- 
tions will be strongly model dependent. 
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